Sub-cycle waveform synthesis for effective control of electron localization 
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We have proposed and demonstrated an effective scheme to steer both the electronic and 
molecular motions of hydrogen molecule and its isotopes. Our scheme lies in the idea of 
sub-cycle waveform shaping (SWS) with synthesizing coherent IR pulses of different wave- 
lengths. The quantum model simulation indicates a remarkable enhancement of the electron 
localization asymmetry compared with the few-cycle pulse. 68% and 89% electrons of the 
dissociating fragmentation can be localized to one nucleus for H2 and H2, respectively. 

Coherent control of quantum dynamics in chemical reactions and photoelectron processes has 
been an exciting goal of chemistry and ultrafast physics. Such a process naturally involves both the 
electron and nuclear motions. On the one hand, femtosecond chemistries have made remarkable 
advance to control the pathway of molecular fragmentation using a optimally-shaped long pulse 
l|-l5|. On the other hand, the advent of carrier-envelope phase(CEP) stabilized few-cycle (~ 5 fs) 
pulse and the attosecond pulse have represented a robust tool to steer the electronic dynamics on 
the attosecond time scale. One interesting question therefore is whether can we combine these two 
directions by effectively steering both the electronic and molecular motions. 

It was demonstrated that, by using a few-cycle pulse, the electron wavepacket can be localized 
to one specific nucleus and therefore shows asymmetric electron localization for hydrogen molecules 



(or its isotopes) [6|-l8(] and the CO 



Another approach to localize the electron was also proposed 
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and demonstrated by synthesizing of fundament laser with another extreme ultraviolet 



infrared pulse 
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131] or 



14l |. However, the localization asymmetry of the total dissociation in the few- 



cycle pulse is quite small [15|] even through the kinetic energy release (KER) spectrum reveal a 
large asymmetry. If the molecular wavepacket could respond to the electronic motion, it would be 
possible to control the electron location and chemical band in a very direct way. Nevertheless, the 
precise steering of electronic motion requires the control of sub-cycle waveform, thus an extreme 
short pulse, e.g., 3.5-fs 



la ], is desired. On the contrary, the vibration period is longer than 14 fs 
even for the simplest molecule ion, i?^. The nuclear motion is too slow to respond directly to the 
electronic motion and the few-cycle laser pulse. To meet this challenge, in this letter, we proposed 
and numerically demonstrated a novel scheme to effectively steer the electronic and molecular 
motions of hydrogen molecule and its isotopes. Our scheme lies in the idea of sub-cycle waveform 
shaping (SWS) with synthesizing coherent IR pulses of different wavelengths. We show that the 
SWS enables us to significantly enhance the asymmetry of electron localization compared with the 
few-cycle pulse. 

With the development of laser technique, it becomes possible to synthesize two coherent light 



sources output from different fiber laser 



17!], optical 



parameter amplifier (OPA) 18 



191 ] or optical 



201 ]. The basic idea of SWS stems from the 



parameter chirp pulse amplifier (OPCPA) systems 
temporal profile modulation by the beat wave of two IR waves 211]. By stabilizing the phase of 
the IR pulse, the waveform of the synthesis pulse can be precisely control. As shown Fig. [l]||a), 
except the weak pedestal peaks, the synthesis field with a 15-fs Ti:saphire laser and 25-fs IR (1200 
nm) laser is equivalent to that of a 4-fs one-color (OC) field, which thus enables to precisely 
control the electronic motion. Moreover, the synthesized pulse contains pedestal peaks maintain 
a duration of about 15 fs, which is comparable to the vibration period of H2 and nuclear motion 
can respond to these pedestal peaks. Therefore, we can expect an effective control of both the 
electronic and molecular motions with this SWS field. We need to emphasize that to precisely 
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FIG. 1: (Color) (a) Temporal profiles of laser amplitude of the 4-fs fundamental (800 nm) OC pulse (red 
dashed line) and SWS pulse (blue line, 15 fs/800 nm + 25 fs/1200 nm). The green filled line shows the 
ionization rate of H2 subjected to the SWS field. The laser intensity is 1 x 10^'^ W/cm^ with an intensity 
ratio of 15%. 0ce and are 0. (b) CEP dependence of the localization asymmetry. The solid magenta and 
blue lines correspond to the interaction of SWS pulse with 4th and 8th vibrational states of , respectively. 
The dashed red and green lines correspond to the interaction of 5-fs OC pulse with 4th and 8th vibrational 
states of H2 , respectively, (c) The CEP dependence of the Franck-Condon averaged localization asymmetry 
for the 5-fs OC field (red dashed line) and SWS pulse (solid blue line), respectively. 



control the electronic dynamics, it is important to shear the waveform of SWS pulse similar to a 
few-cycle OC pulse. This condition requires that the pulse duration of the synthesized IR pulse 
should be less than 25 fs. 

To demonstrate our scheme, we use a numerical model that has been well established for study- 



ing electron localization with few-cycle pulse 



7,E 



10 



12 



14l | . This model solves the time-dependent 



Schrodinger equation for the nuclear wavepacket in the Born-Oppenheimer (BO) representation. 
The time scale of the rotation is several hundreds fs, thus can be neglected in the ultrashort laser 
pulse. We concentrate on the hydrogen molecule and its isotopes, which has been of special in- 
terest for extensive investigation. The molecular dynamics can be well described in terms of the 



two lowest-lying electronic state. There are generally two kinds of experiments, one concentrates 



on the target of H2 a, 



23 



12l |. The underlying physics is 



24l | and the other one on the H2 7, 
however very different. We first analysis the first situation and assume that H2 is aligned along 
the direction of the laser field. Then the TDSE can be written as 7| 



(2m) + 



gu 



+ vj UJ' 



9 9U 



(1) 



where /i is the nuclear reduced mass, Vg and Vu are the BO potential curve calculated from the 
electronic structure of the lag and 2pau states of H2 and the coupling Vgu = E{t)Dgu with Dgu 
being the the electronic transition dipole between these two states. The electric field is expressed 
by ^mix = ^oexp[-2 In 2(^)2] cos(woi + 0ce) + ^1 exp[-2 In 2(^)2] cos{Kujot + (/)ce + 0i), where 
K is the frequency ratio between the main fundamental field (wq) and the weak mixing IR field 
(wi = KujQ, K < 1), the subscript (0, 1) means the two laser field components, E is the electric 
field, r, (j)cE and 0i denote the pulse duration (FWHM), the CEP and relative phase, respectively. 
Here (pi = and the CEPs of the fundamental and the IR signal lights are equal due to the phase 
relationship of OPA ^]. The probability of electron localized on the left or right nucleus can 
be calculated by introducing the left and right localized wavefunction ipi^r = {^g i ^u)\/2- The 
kinetic energy release (KER) spectrum Si^r{Ek) is calculated by Fourirer analysis of the real-space 



wavefunction with the bound states being projected out 



25l |. The electron localization is quantified 



by the left-right ejection asymmetry parameter A{Ek,(l)cEp) = gjt'.feEpi+Sjt-i^cEp) ' 

Figure HJb) shows the electron localization asymmetry as a function of CEP (pcE- The solid 
magenta and blue lines correspond to the interaction of SWS pulse with 4th and 8th vibrational 
states of H2, respectively. For comparison, a similar simulation has been performed for the 5-fs 
OC field and the results are presented by the dashed red and green lines in Fig. d^b). In both 
situations, we can observe sensitive CEP dependence of the electron localization. Nevertheless, in 
comparison with the 5-fs OC pulse, the asymmetry amplitude in the SWS pulse is two- and three- 



fold higher for the 8th and 4th vibrational states, respectively. Moreover, we need to emphasize that 
preparation mechanism of H2 generally indicates an incoherent Franck-Condon(FC) distribution 



of vibrational states in experiments 



231. 



vibrational states weighted by FC factors 



'hus we have to average the observables over the initial 



251]. Figure [D^c) shows the FC averaged localization 



asymmetry. However, we can see very distinct situations: in the OC field, the asymmetry has 
been significantly smeared. By contrast, very large asymmetry with amplitude of 0.35 still can be 
obtained in the SWS pulse, which corresponds to an electron localization probability about 68% 
of all fragmentation events. We have changed the intensity from 2 x 10^^ to 2 x 10^^ W/cm^ and 
pulse duration from 3 to 10 fs, the FC averaged asymmetry in the OC few-cycle pulse is below 
0.06, which is only 1 sixth of that in the SWS pulse. 




FIG. 2: (Color) KER spectra as a function of CEP in the (a) SWS pulse (b) 5-fs pulse. Asymmetry map 
as a function of KER and CEP in the (c) SWS pulse and (d) 5-fs pulse. Asymmetry map as a function of 
vibrational state and CEP in the (e) SWS pulse and (f) 5-fs pulse. The laser parameters are the same with 
Fig. mh). 
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Figure [5] shows the distinct difference of the electron locahzation between the SWS [(a), (c) 
and (e)] and OC 5-fs pulses [(b), (d) and (f)], respectively. From these results we can get several 
important remarks, (i) KER is lower in the SWS pulse and concentrated around 0.5 eV. The 
higher KER (> 1 eV) is significantly suppressed compared with that of OC 5-fs pulse [see Fig. 
[2ja) and (b)]. (ii) The asymmetry electron localization depends sensitively on the CEP and KER 
in the OC 5-fs pulse, therefore resulting a "tilted stripe" as shown in Fig. ^d). However, in the 
SWS field, the electron localization depends only sensitively on the CEP, but not sensitively on 
the KER. Therefore, as shown in Fig. EJc), we observe a "vertical stripe". Note that the total 
population of the electron localization is the integration over the KER. Consequently, OC 5-fs 
totally results in a very weak contrast of the electron localization. Nevertheless SWS pulse can 
maintain the large contrast, resulting a much more effectively control of electron localization, (iii) 
The CEP dependence of electron localization is "in-consistent-step" for different vibrational states 
in the SWS field, however that is "out-of-step" in the OC field. As shown in Fig. [2l|e), a "tilted 
stripe" can be observed in the asymmetry map as a function of CEP and vibrational states. Thus 
electron localization of some vibrational states (e.g., v=4 and 8 shown in Fig. [Hb)) are completely 
reverse, therefore smearing out the electron localization. However, in the SWS pulse, more electron 
is localized in the left nucleus in a broad range of CEP (0.27r to O.Tvr) for the lowest 10 states. 
Note that the higher vibrational states correspond to a EC factor 1 order of magnitude lower, thus 
play a minor role. Such results clearly demonstrate the effective control of the electron localization 
with the SWS pulse. The FC-averaged asymmetry parameter is one order of magnitude higher 
than that of OC 5-fs scheme. Moreover, we should not that to keep the waveform of SWS pulse 
similar to a few-cycle pulse (< 5/s), the pulse duration of the IR pulse should be less than 25 
fs. Otherwise, different KER will show a different CEP dependence but the total electron 
localization asymmetry will be significantly washed out. 

In order to shed a deeper insight onto the difference of the electron localization in the SWS 
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FIG. 3: Color)Top plot: evolution of the eigenenergy of the quasi-states (green line) and (/)2(blue line) 
and the intcrnuclcar distance (red line). Bottom plot: evolution of the left (red line) and right (blue line) 
localized wavefunction | (j)i^r (a), (b) and (c) correspond to the interaction of SWS pulse with the state 
v=3, 4 and 8, respectively, (d) and (e) correspond to the interaction of 5-fs pulse with the state v=3 and 8, 
respectively. The CEP (t>ccp = OAtt and other parameters are the same with Fig. [TJ 



and OC few-cycle pulses schemes, we have analysis the molecular dynamics in terms of the quasi- 
static states, which is related to left and right localized wave functions 'i/'i.r via 



261, 



27[, -01,2 



[-0; (cos 6* lb sin6') it il^ricosO =F sin^)], where 29 = tan ^{—2Vgu/{Vu — Vg)). The corresponding 
quasi-static eigenvalues are Vi,2 = -^^^ii^ =F ^J{Vu-Vg)'^/4 + V^^. For small internuclear distance 
R < 3 au, the energy gap between Vi and V2 is much larger. Thus the coupling between (j)i^2 is 
difficult, electron dynamics is adiabatic and can not be localized. As R increase larger than 3.5 au, 
significant nonadiabatic coupling becomes possible and electron localization will be established. If 
R > 6 au, the internuclear barrier is very high, which frozen the electron localization. Figure [3] 
presents the evolution of the quasi-state eigenvalues Vi^2 and internuclear distance R (top plots) and 
the left and right localization probability (bottom plots) for different vibrational states subjected 
to the SWS and 5-fs OC pulses, respectively. Here the CEP cpcE is set to 0.4 tt [the corresponding 
asymmetry map is marked by the arrow in Fig. and (f)]. For the interaction of lower vibrational 
state with 5-fs OC pulse, as shown in Fig. [3jd)(v=3), the motion of nuclear wavepacket is too slow 



to respect the ultrashort laser pulse. Thus the internuclear distance slightly changes around the 
equilibrium 2 au. The electron follows the oscillation of the laser field and changes its localization 
in each half cycle. The dynamics for higher vibrational states is completely different. As shown 
in Fig. [3]^e)(v=8), the energy gap becomes smaller. The internuclear distance gradually increases 
larger than 6 au at the end of the pulse and more electron is localized near the left nucleus. By 
contrast, in the SWS pulse, the pedestal pulse is effective to stretch the internuclear distance, which 
increases to 4.5 au for V=3 and larger than 6 au for the higher states. Therefore, these vibration 
states show a similar nonadibiatic behavior. Of course, the CEP dependence can not be exactly 
same for different states, the different nevertheless becomes much smaller compared with that in 
the 5-fs pulse. This is possibly the reason why the CEP dependence is "in consistent step" for the 
SWS pulse. Moreover, because the internuclear distance is stretched, the dissociation happens at 
a large internuclear distance. Consequently, the KER is lower and the high KEP is significantly 
suppressed compared with that of the 5-fs pulse. 

Another advantage of our scheme of synthesis the coherent IR laser pulses output from OPA 



Ig or OPCPA systems 



201 1 is that the central wavelength of the ideal or signal lights can be 
adjusted in a broad range. This provides additional freedoms to control the sub-cycle waveform 
and the electron localization. Figure SH^a) shows the CEP dependence of asymmetry of electron 
localization by changing the IR wavelength to 1300 nm (blue dashed line), 1500 nm (red solid line) 
and 1800 nm (green dash-dotted line), respectively. A large asymmetry amplitude of about 0.3 
can be obtained in all these cases, which corresponds to 65% localization probability of the total 
dissociation events. We have also simulated the electron localization for D^. The CEP dependence 
of localization asymmetry is shown in Fig. Hl^b) for the SWS field synthesized with 800-nm and 
1200-nm (blue dashed line), 1500-nm (red solid line) IR pulses, respectively. The laser intensity is 
1 X 10^"^ W/cm^, relative intensity is 15% and all the pulse durations are 25 fs. The localization 
asymmetry shown in Fig. Hl^b) demonstrates the effectiveness of our SWS scheme in the case of 
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FIG. 4: (Color) (a) CEP dependence of the asymmetry of electron localization of H2 subjected to the pulse 
synthesized by 800-nm and 1300-nm (dashed blue line), 1500-nm (solid red line) and 1800-nm (dash-dotted 
green line), respectively. The pulse duration of the 800-nm pulse is 15 fs and is 25 fs for the other IR pulses. 
The laser intensity is the same with Fig. [TJ (b) Same to (a), but the target is changed to 13 J and the pulse 
duration of the 800-nm pulse is 20 fs and it 25 fs for the other IR pulse, (c) Same to (a), but the target is 
changed to H2- The pulse duration of the 800-nm pulse is 15 fs and is 25 fs for the other IR pulses. The 
laser intensity is the same with Fig. [TJ 



3, 191]. Briefly, this 



A different mechanism of electron localization lies on the recollision process 
process starts from the ionization of deuteron molecules. The electron is accelerated in the laser 
field and recollides with the molecule ion as the electric field changes its direction. The recollision 
can excite the electron wavepacket and may localize the electron to one nucleus. Different from 
the previous process, where the electron localization is induced by the laser coupling between two 
states, the recollision process plays the vital role for electron localization. The SWS pulse also 
enables us to effectively control the recollision and so the electron localization. To demonstrate our 
scheme, we consider the target of i?2- As we first we calculated the ionization of H2 with the 
molecular Ammosov-Delone-Krainov (MOADK) ionization model [28]. The green filled area in Fig. 
[TJa) shows the ionization rate of H2 subjected to the SWS pulse. One can see that the waveform 
of the SWS is similar to a 4-fs pulse. The electron mainly ionized at the highest peak, while the 
ionization rates at other peaks are less than 5% of the highest peak. At each ionization time, the 
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electron wavepacket of H2 is launched to the lscr„ state of H2 ■ This wavepacket propagates for 
0.7To, then the electron recollides with the H2 and boost the wavepacket to 2pau state, which 
is repulsive and can be quickly dissociated. We use the same two-state model of TDSE to follow 
the motion of the wavefunction. Finally the resulting localization probability summed incoherently 
with the weight given by the ionization rate over all the wavepacket launched at different ionization 
times. Figure m^c) shows the asymmetry of electron localization in the SWS pulse synthesized by 
800-nm and 1200-nm (blue dashed line), 1500-nm (red solid line) IR pulses, respectively. One can 
see a very large asymmetry of electron localization, especially in the later case, the asymmetry 
amplitude is 0.78, which corresponds to an electron localization probability as high as 89% of all 
fragmentation events. We should emphasize that the waveform of our SWS pulse is similar to a 
4-fs pulse, which results in solely one ionization and recollision events. Such a character is crucial 
to effective control of electron localization, because the asymmetry maybe washed out if multi- 
recollision events contribute [7|. Note also that Ref. [lO(] recently reported a relevant experiment, 
which demonstrate that, for molecule the electron localization induced by laser coupling also 
can be dominant. This process starts from the ionization of H2 and then the wavepacket is coupled 
by the laser field but the recollision process plays a minor role. We have also simulated this 
process with the SWS pulse. Our results have indicates a high degree of control over the electron 
localization (asymmetry parameter is about 0.4). 

In conclusion, we have shown an effective scheme based on the SWS by synthesizing multicycle 
IR pulse of different wavelengthes to control the electron localization of H2 , H2 and its isotopes. 
Compared with the few-cycle pulse, SWS pulse shows a higher degree of electron localization 
control. In the case of H2 , a localization probability of 68% are realized and for H2, the localization 
possibility can be as high as 89% with the SWS scheme. Moreover, the current laser technology has 
made it possible to synthesize and stabilize the phase of multicolor coherent pulses, e.g., the CEP- 
stabilized pump, signal, idle lights of optical parametric synthesizers. We may expect an effective 
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control of the electron localization for more complex molecules than hydrogen by multicolor SWS. 
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